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ABSTRACT 

The experimental stress intensity factors for various chevron 
notched four point bend specimens are presented. The experimental 
compliance is verified using the analytical solution for a straight 
through crack four point bend specimen and the boundary integral 
equation method for one chevron geometry. Excellent agreement is 
obtained between the experimental and analytical results. In this 
report, stress intensity factors, loading displacements and crack 
mouth opening displacements are reported for different crack 
lengths and different chevron geometries, under four point bend 
loading condition. 


INTRODUCTION 

The usefulness of the chevron notched specimens for testing 
the fracture toughness of brittle material is well established 
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(Ref. [l]-[4]). However, the chevron notch bend specimen has not 
been analyzed as rigorously as the chevron notched short rod and 
bar specimens, ([5]-[12]). Munz et al. [13] used Bluhm's slice 
model [9] to obtain the stress intensity factors (SIF) for the four 
point bend specimen while Wu [14] used the same method to obtain 
the SIF for the three point bend loading condition. Shih [15] 
attempted to match a valid fracture toughness K, c value of a 
specific material with the maximum load to failure of a three point 
bend chevron notched specimen. To the author's knowledge, the only 
rigorous analysis for the three point bend specimen has been 
conducted using the boundary integral method [16] and the finite 
element method [17], Furthermore the only experimental compliance 
calibration of the bend specimen has been the original work of 
Bluhm [9], [18] and Oucherlony [19] for the three point bend round 
bar. The continued interest in the chevron notched four point bend 
specimen especially for elevated temperature testing has prompted 
the current derivation of the stress intensity factor from 
experimental compliances for various chevron notched four point 
bend specimens . 

COMPLIANCE MEASUREMENTS 

Experimental Procedure 

The four point bend loading fixture is shown in Fig. 1. This 
fixture allows the bottom support rollers and the top loading 
rollers to roll freely to eliminate errors due to any sliding 
friction. The load is transmitted through a conical pad on the top 
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fixture to minimize any misalignment in the rollers and/or non- 
parallelism in the specimen. The upper fixture is lifted off and 
lowered onto the specimen through a spring system that is detached 
when the compliance measurements are taken, Fig. 1. The load 
displacement is taken as the relative vertical displacement between 
small pins located along the centerline of the specimen directly 
above the support rollers and below the loading rollers. Fig. 2-a. 
The motion is measured with an LVDT by attaching the core of the 
LVDT to a "T" shaped hanger that rests on the pins below the 
loading rollers, while the LVDT body is secured to a saddle 
arrangement resting on the pins above the support rollers. The 
support span ( S 1 ) is fixed at 8" while a loading span (S 2 ) of either 
2" or 4" is used. The roller diameters are changed to match the 
specimen according to ASTM E399 recommendations [21] for three 
point bend tests. The specimens are machined from a two inch thick 
plate of 7075-T651 aluminum in the L-T orientation (ASTM Coding) to 
a thickness B, of 1”, a width W of 1" or 2”, and a length of 8.25". 

The chevron notches have a relative notch length at the 
surface (a^W^,) of 1.0, and a relative notch length to the chevron 
tip (a 0 /W=a 0 ) of 0.0, 0.2, and 0.5 (see Fig. 2-b) . The width of the 
chevron notch is 0.035" with a semicircular bottom. The notch is 
incrementally extended from a 0 to an a (a=a/W) of 0.95 by sawing a 
0.022" notch on-a thin blade band saw. Load-deflection traces are 
made after each incremental extension of the notch. The load is 
progressively decreased as the notch depth is increased to keep the 
applied stress intensity factor below lOksiyin. 
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Analytical Procedure 

To support the experimental compliance method, a three 
dimensional elastic analysis is performed on one chevron notched 
geometry using the boundary integral equation method (BEM) . This 
method uses only boundary surface elements. The BEST3D computer 
code [22], which employs a quadratic variation of the boundary 
traction and displacement over triangular and rectangular elements, 
is used. Due to the symmetries in the loading and the geometry, 
only a quarter of the specimen was modelled. The quarter of the 
beam is also divided into three subregions that are interconnected 
by interface elements since the aspect ratio of the beam length 
over the beam width is large, Fig. 2-c. Making use of the half 
symmetry option in the BEST3D code, surface elements were not 
required in the YZ plane at x=0, as seen in Fig. 2(c). The shaded 
area along the XY plane represents the un-cracked ligament of the 
chevron notch for a=a/W=0.5. The finite width slot N cut into the 
specimens to form the chevron is not modelled. The load line 
displacements are calculated for one chevron geometry (with 
a 1 =a 1 /W=l . , a 0 =a(j/W=O . 2 and B/W=0.5) and different crack length a/W 
ranging from 0.25 to 0.75 for two span length ratios (S 1 /S 2 =2. and 
4 . ) . Crack mouth opening displacement is calculated for different 
crack increments and compared to the experimental results . The 
analytical calibration of the stress intensity factor K simulated 
the experimental calculations in determining the change in the 
external work as the crack advances: 
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( 1 ) 





Where G is the energy release rate 


E is the elastic modulus 


P is the total applied load 


D is the load line deflection 


A is the cracked area 


Before performing the experimental calibration on the chevron 
notched four point bend specimens, the accuracy of the experimental 
set up is tested by calibrating a straight through four point bend 
specimen and comparing it with published results. 

RESULTS 

Straight Through Notch Specimens 

The experimental load line compliance for the straight through 
notched four-point bend specimen is normalized in the following 
manner, to cancel the effects of various span lengths: 


where B is the beam thickness, W is the beam height, and C 0 is the 
un-notched non-dimensional beam compliance obtained from Bluhm [9]: 



( 2 ) 
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(3) 
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4 W 2 (S 1 - S 2 ) 


where u is the Poisson's ratio. 

The normalized load line compliance D N is plotted as a 
function of the normalized crack length a in Fig. 3 for different 
span ratios and specimen geometries. The normalized compliance D N 
is fitted using the least squares method to give: 


„ _ 9 .305510 + 4.30643a 2 - 4.55453a 3 - 1.02783a 4 + 2.04204a 5 , x 

d n — Tj (4) 

(1 - a) 2 

The experimental data as well as the values from Eq. 4 are 
shown in Table 1. A good fit is observed over the full range of a= 
0 to 0.95, independent of span lengths and beam heights, as seen in 
Fig. 3. 

The normalized stress intensity factor Y is determined based 
on Eq. 1 as the change in specimen compliance with increased crack 
length: 


K_ B t? 3/2 
P < S 1 - S 2 ) 


'I 5 
2 6a 



(5) 


The calculated normalized stress intensity factors Y are 
listed in Table l with those from the Srawley-Gross relationship 
for pure bending [23] given here for completeness: 
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Y = — 


2 \ 


(l-O ) 3 


( 6 ) 


1.9887 -1.326a - 


(3 .49 -0 ,68a +1 ,35a 2 ) a (1-a) 
( 1 +a ) 2 


The agreement between the experimentally derived normalized 
stress intensity factor Y and the analytical solution is quite good 
over the full range from a=0 to 0.95 as shown in Fig. 4. As Eq. 6 
indicates, the normalized stress intensity factor is independent of 
the specimen geometry and the span ratios. 

Chevron Notch specimen 

The normalized load line compliances D N for the chevron 
notched four-point bend specimens are also calculated using eq. 2 
which was shown to render the D N values independent on the span 
lengths and the beam height used. The un-notched dimensionless 
compliance is again computed from Eq. 3. The normalized load line 
displacements are shown in Fig. 5 along with the analytic solution 
from the BEM for a chevron notch specimen with an a=0.2 and a^— 1.0. 
The values from the analytical solution are somewhat lower than the 
experimental values. This is due in part to the difference in the 
value used for the un-notched compliance and also since the 
analytic solution is based on a sharp crack while the actual 
specimens had a finite width slot. The method of least squares was 
used to fit an expression to the normalized load line compliance 
for each notch geometry. These expressions had the same form as 
Eq. 2 : 
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Ar 



(1-a) 2 


where the least square coefficients f 5 are listed in Table 2 for 
different crack geometries. 

The normalized stress intensity factor Y* for a chevron notch 
is obtained from the following relation given by Munz [13]: 


Y* 


KB W 3/2 


1 

a i~ a o 6 ' 

P(S 1 -S 2 ) 


2 

a-a 0 5a N 


( 8 ) 


The values of Y* are plotted in Fig. 6 along with values from 
the analytical solution for a 0 =0.2, a.,=1.0. Good agreement between 
the values derived from the data and the BEM results is observed. 

Crack Mouth Opening Displacement 

In addition to the load line deflections, the crack mouth 
opening displacement V is determined with an ASTM E399 type clip 
gage. The crack mouth opening displacement (CMOD) measurements are 
simpler and more direct then the load line displacement, thus they 
lend themselves very well to automated cyclic crack growth and "R" 
curve studies. The normalized form of the crack opening 
displacement is formulated in the following manner: 



where V N is now normalized with respect to the moment P(S 1 -S 2 ) and 
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modulus , E . 

The experimental normalized CMOD data as well as the 
analytical values is plotted in Fig. 7. The values from the 
analytical solution are somewhat lower. This is due in part to the 
difference in the notch geometry. In the analytical solution the 
notch width is zero while the test specimens have a notch width of 
0.035". The method of least squares is also used to fit an 
expression of the same form as Eq. 3 to the normalized CMOD data: 


EVB _ ^1 

( . w ) 

f 0 + f 1 a+f 2 a 2 + f 3 a 3 + f 4 a 4 + f 5 a 5 

P 0 

Ul -Sj 

(1-a) 2 


where the least square coefficients f { are given in Table 2. 

In crack growth and R curve studies, determining the crack 
length from the crack opening displacement is more convenient. 
Therefore, inverse expressions were obtained following the form 
used by Haggag and Underwood [24], to estimate the crack length 
given the CMOD value: 


a = f 0 + f x U + f 2 U 2 + f 3 U 2 + f 4 U 4 + f 5 U 5 (11) 


where 


U = 


1 

* i 


The least square coefficients f { for different crack geometries are 
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also given in Table 2. 


CONCLUSIONS 


1. The excellent agreement between the Srawley-Gross analytical 
results and the experimentally derived stress intensity factors for 
a straight through crack in a four point bend specimen 
substantiates the method used in obtaining the experimental load 
line displacements. 

2. The experimental data confirms that the normalized stress 
intensity factor Y is independent of span ratio and specimen 
geometry for the straight through crack four point bend specimen. 

3. The experimental study performed in this study provides the 
necessary numerical results to determine the stress intensity 
factor for several four point bend chevron notched geometries. 
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Table 1- Normalized load line displacements and stress intensity 
factors for four point bend straight through cracks (E=10.5MSi, 
l )=0 . 3 ) 


W,in 

1.002 

1.984 


B,in 

1.001 

1.001 

S r in 

8.00 

8.00 

S 7 ,in 

2 

4 
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a 


a 

0 N 

a 

Y 


DATA 

EQ.4 


DATA 

EQ.4 


EQ.5 

REF [ 23 ] 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 

- 

0.093 

0.056 

0.073 

0.076 

0.099 

0.081 

0.084 

0.05 

0.65 

0.64 

0.195 

0.31 

0.21 

0.29 

0.203 

0.41 

0.32 

0.10 

0.85 

0.89 

0.298 

0.72 

0.70 

0.71 

0.298 

0.66 

0.71 

0.15 

1.05 

1.08 

0.352 

1.08 

1.07 

1.04 

0.347 

1.00 

1.01 

0.20 

1.24 

1.26 

0.400 

1.43 

1.33 

1.43 

0.398 

1.35 

1.41 

0.25 

1.43 

1.44 

0.451 

1.99 

1.95 

1.96 

0.457 

2.06 

2.04 

0.30 

1.64 

1.64 

0.504 

2.74 

2.67 

2.70 

0.505 

2.85 

2.72 

0.35 

1.87 

1.86 

0.554 

3.60 

3.57 

3.65 

0.553 

3.61 

3.63 

0.40 

2.13 

2.12 

0.602 

4.94 

4.78 

4.92 

0.604 

4.95 

4.97 

0.45 

2.44 

2.43 

0.656 

6.88 

6.91 

6.99 

0.651 

6.55 

6.75 

0.50 

2.81 

2.81 

0.698 

9.54 

9.32 

9.42 

0.695 

9.67 

9.19 

0.55 

3.28 

3.30 

0.750 

14.50 

14.25 

14.26 

0.753 

14.58 

14.61 

0.60 

3.90 

3.94 

0.800 

22.87 

22.52 

22.97 

0.806 

24.39 

24.44 

0.65 

4.74 

4.82 

0.859 

47.86 

47.37 

47.44 

0.856 

44.68 

45.67 

0.70 

5.96 

6.06 

0.901 

100.6 

98.29 

98.86 

0.902 

99.62 

101 .6 

0.75 

7.84 

7.95 

0.946 

364.2 

359.0 

350.7 

0.955 

470.1 

503.5 

0.80 

11.01 

11.10 








0.85 

17.09 

17.06 








0.90 

31 .76 

31.33 








0.95 

91.15 

88.69 
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Table 2- Least square fit for the experimental calibration of chevron notched four point 
bend specimens (E=10.5MSi, i>=0.3) 


NORMALIZED LOAD LINE DISPLACEMENT 

D 

-f EDB _ 

Co 

{ W V f 0 + / 1 a + f 2 a 2 +f 3 a 3 +f 4 a 4 +f 5 a 5 

N [ p 

[S 1 ~ S 2 ) (1-a) 2 

°n 

a 1 

W/B 

_ fn ... 

fi 


f* 

. _ f 4 

f K 

0 

1 

2 

0.4157 

-0.6514 

6.2589 
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1.8503 
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0.2 

1 
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1.3099 

-5.8070 
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-28.178 

17.932 

-4.1748 

0.5 

1 

2 

0.8610 
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349.35 

-302.48 
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0 

1 

1 

0.554 

-1.4107 

8.8164 

-14.661 

10.412 

-2.6803 

NORMALIZED CRACK MOUTH OPENING DISPLACEMENT 

II 

’ EVB ' 

( W f 0 +f 1 a+f 2 a 2 + f 3 a 3 +f 4 a 4 + .f 5 a 5 

P 

\ s i ~ s 2 j (1-a) 2 

g n 

“1 
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f n 

fi 
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^A 

f c 

0 

1 

2 

1.6958 

-0.2878 

16.144 

-33.834 

29.0289 

-8.7674 

0.2 

1 

2 

4.7670 

-18.477 

62.901 

-96.074 

71.914 

-21.110 

0.5 

1 

2 

121.922 

-754.42 

1927.9 

-2452.4 

1546.83 

-385.633 

0 

1 

1 

2.0193 

-2.8288 

26.5249 

-55.0951 

48.3189 

-14.7658 

NORMALIZED CRACK LENGTH AS A FUNCTION OF THE NORMALIZED CMOD 

a = f^f 2 u+f 2 u 2 +f 3 t7 3 + f 4 £7 4 +f 5 tr 5 
where U = = 

* i 
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W/B 

fn 
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f ? 

f. 
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fi 

0 

1 
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1.01069 
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0.2 

1 

2 

1.01269 

-2.8250 
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-140.298 
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-683.339 

0.5 

> 1 
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-0.54499 

-30.981 
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-576.58 

-328.473 

0 

1 

1 

1.02438 

-3.3387 

21.2819 

-147.101 

417.057 

-419.430 
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Table 4- Normalized stress intensity factor Y for chevron notched four point bend specimens (E=1 0.5MSi, 
u=0.3) 
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9.55 
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.85 
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20.33 

18.64 
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33.10 

33.60 


35.23 

33.47 
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92.19 


96.51 
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Table 5- Normalized crack mouth opening displacement for chevron notched four point bend specimens (E=10.5MSi, u=0.3) 
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Figure 1.— Four point bend loading fixture. 
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Section view 


Side view 



(b) 


Figure 2. — (a) Schematic of a four point bend chevron notch specimen, (b) Detail view of the chevron notch geometry. 
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Figure 2. — Concluded, (c) Typical boundary integral mesh with quarter symmetry. 



0 0.2 0.4 0.6 0.8 


Normalized crack length a/W 

Figure 3. — Normalized load line displacements for a straight 
through crack as a function of the crack length. 
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Figure 4. — Normalized stress intensity factors for a straight 
through crack as a function of the crack length. 



Figure 5. — Normalized load line displacements for various 
chevron notch geometries as a function of the crack length. 
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Figure 6. — Normalized stress intensity factors for various 
chevron notch geometries as a function of the crack length. 



Figure 7. — Normalized crack mouth opening displacements for 
various chevron notch geometries as a function of the crack 
length. 
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